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SUMMARY 


A  consistent  estimate  is  proposed  for  the  scale  parameter  /f  in 

n.KyH. ’ ‘  e s'-  -  mu 

the  model  Y.  =  y.  +  e  ,  Is  i  £  n,  where  the  are  unknown  location 

iu.1"  1  1  1  v,.  /  p  »«•  1 

c 

parameters  and  the  e^  are  independent,  identically  distributed  random 
errors  with  density  function  f.  This  parameter  arises  in  the 
variance  formula  for  rank  estimates  of  location.  The  proposed  estimate 

is  based  on  differences  of  residuals  Y- -y  • ,  where  y  •  is  an  estimate 

'  -  h  ■ 

of  ,y^.  When  the  yu  follow  the  structure  of  the  general  linear  model, 
the  estimate  is  shown  to  be  consistent  under  the  usual  assumptions  on 
the  design  matrix.  The  estimate  does  not  require  the  symmetry  of 
the  density  f. 
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1.  INTRODUCTION 


Let  n  observations  Yj^ . Yn  follow  the  model  , 

Is  is  n,  where  the  y^  are  unknown  location  parameters  and  the  e^  are 
independent,  identically  distributed  random  variables  with  density  f. 

Of  particular  interest  here  is  the  case  where  the  y^  follow  the 
structure  of  a  linear  model. 

Consider  the  problem  of  estimating  the  scale  parameter 
y  *  /f  (x)dx.  This  parameter  arises  in  various  statistical  procedures  based 
on  ranks.  For  instance,  in  one -sample,  two-sanple  and  linear  model  problems 
it  appears  in  the  variances  of  rank  estimates  when  Wilcoxon  scores  are  used. 

In  linear  model  problems,  Hettmansperger  and  McKean  (1976),  (1977)  proposed 
test  statistics  that  are  based  on  the  drop  in  a  dispersion  function  between 
full  and  reduced  models  and  this  difference  must  be  scaled  by  y. 

For  the  one-sample  and  two-sample  problems,  Lehmann  (1963)  proposed  a 
consistent  estimate  of  y  based  on  the  length  of  a  di s t r ibut i on - free  confidence 
interval  for  a  shift  parameter.  For  the  linear  model  problem,  under  the 
assumption  that  f  is  symmetric  about  zero,  Hettmansperger  and  McKean  (1976), 
(1977)  proposed  using  Lehmann's  estimate  of  y  computed  from  the  residuals 
and  proved  it  was  consistent. 

In  this  paper  an  estimate  of  y  is  proposed  that  is  applicable  for  linear 
models  and  more  general  situations  without  assuming  symmetry  for  the  density  f . 
To  describe  this  estimate  consider  residuals  -  y^,  where  is  an 

estimate  of  y^,  Is  is  n.  Let  Hn(t)  be  a  weighted  empirical  cdf  of  the 
absolute  values  of  differences  of  residuals, 


I 
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H At)  =  l  .  a. . 4>( I Z- -Z.  | ,  t) 
n  l<j  J-J  J  1 


where  4>(u,v)  =  0,  *5,  1  as  u  >  v,  u  =  v,  u  <  v  and  a^  >  0  with 


l  ■  »• 

1<J  J 


The  proposed  estimate  of  y  is 


(1.1)  y  =  fin(t/»€  )/(2t/tfi  ) 

th  A 

where  t  is  the  a —  quantile  of  Hn(t)  ,  0  <  a  <  1. 

Another  similar  estimate  is 

(1.2)  y*  =  $n(Wk)/2Wk  =  k/2MVk 

where  s  . . .  <  are  the  M  =  (”)  ordered  absolute  differences 

J  Z . -Z - 1  for  i  <  j  and  k  is  chosen  so  that  W,  *  0  (1) . 

'Jl  K  P 

The  unweighted  (or  equal  weight)  case  provides  the  basic  estimate 

while  weights  can  be  used  for  special  purposes.  For  example ,  the  use 

of  0-1  weights  would  restrict  the  number  of  differences  used  in  H  and 

n 

reduce  computational  effort.  In  analysis  of  variance  models  and  stratified 
sampling  models  where  the  data  arise  naturally  in  groups,  the  use  of  weights 
a^j  =  0  if  Zj,  2^  are  from  different  groups  would  leave  the  estimate 
dependent  only  on  the  within -group  variation  of  the  data.  This  intuitively 
appealing  estimate  has  apparently  not  been  considered  in  the  literature. 

For  the  one-sample  problem  with  y^  »  y  for  all  i  and  under  the  assumption 
of  a  symmetric  distribution,  Antille  (1974)  considered  an  estimate  of  y  based 
on  the  number  of  averages  (Zj+Z^/2  in  the  interval  (0,  1/Jn  )  for  i  <  j. 

This  estimate  is  basically  a  window  estimate  of  the  density  of  (¥^+¥2)72 
at  y.  The  window  is  not  symmetric  about  zero  and  the  window  width  does  not 
adjust  to  the  spread  of  the  data. 
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Schweder  (1975)  proposed  a  Lehmann-type  estiruite  of  y  (in  formula 
(3.1))  which  is  basically  a  window  estimate  and  recomnended  a  uniform  window 
symmetric  about  zero.  Assuming  =  0  for  all  i  so  that  no  nuisance  para¬ 
meters  need  be  estimated,  his  estimate  is 

Si,j«iyYii>  v2>/n2hn 

-  C(n-l)/n)Hn(hn/2)/hn  .  (l/nhn) 

i  VY2>'\  > 

where  is  the  width  of  a  uniform  window  and  is  the  empirical  cdf  of 

lY.-Y^]  for  i  <  j.  The  estimate  (1.1)  is  seen  to  be  similar  to  this  type  of 

_w 

estimate  with  a  window  width  of  0  (n  ) .  In  proving  consistency  and  asymptotic 

r 

normality,  Schweder  assumed  the  window  width  satisfies  hn  00  ,  ruling 

-1/3 

out  the  0(n  )  cases,  and  recomnended  hn  =  0(n  ). 

The  estimate  y  differs  from  Schweder 's  estimate  by  allowing  for  the 
estimation  of  nuisance  parameters  in  using  residuals  rather  than  iid 
variables.  Also  it  uses  a  window  width4 estimated  from  the  data.  In  an 
actual  application  of  Schweder 's  estimate  the  choice  of  an  appropriate  s  for 
h^  =  sn*1^3  is  difficult  to  make  and  his  Monte  Carlo  study  indicates  that  the 
bias  of  the  estimate  is  sensitive  to  the  choice  made. 

The  estimate  y*  of  (1.2)  can  be  viewed  as  a  nearest  neighbor  type 
estimate  based  on  the  absolute  values  of  the  differences  of  residuals. 

Nearest  neighbor  density  function  estimates  are  discussed  in  Moore 
and  Yackel  (1976) ,  (1977) . 

Cheng  and  Serfling  (1981)  considered  estimates  of  efficiency-related 


functionals,  including  y  as  a  special  case. 


Hie  main  results  are  presented  in  this  section.  The  proofs  can  be 


found  in  the  Appendix. 

For  each  integer  n  =  1,2,  ...  ,  assume 
(2.D  X  -  ^  +  Sn  . 

where  l n  =  (Yin . Yj n)'  is  an  n  x  1  random  vector,  =  (uln,  ...,  h^)' 

is  an  n  x  1  parameter  vector  and  en  is  an  n  x  1  vector  of  independent, 
identically  distributed  random  variables,  each  with  density  function  f. 

At  times  the  n  will  be  suppressed  to  simplify  the  notation.  For  an  n  x  1 
vector  q  ,  define  an  n  x  1  vector 

2(q)  -  y  -  a 

*  (Z1(a),  ...,  znCD))'  . 

Let  an  indicator  function  be  defined  by  $(u,v)  =  0,*5,1  as  u  >  v, 
u  =  v,  u  <  v.  Suppose  for  each  integer  n  ■  1,2,  ...  a  set  of  weights  is 
given  *  a^  for  1  n  <  j  <=  n  with  a^  =  1.  Consider  the  weighted 

empirical  cdf  of  the  absolute  values  of  the  differences  in  Z(tj)  given  by 

yt.rj.Y)  -  J4<j  . 

This  is  not  quite  a  proper  cdf  since  <J>  ■  \  in  case  of  ties. 

REMARK  2.1 

(i)  1(h)) 

CU)  »n(t.B.Y)lJ4i  H„(t.a-M.Plo 

(iii)  Hn(t,jj(Y),Y)lKi, 

if  jiQf)  is  an  estimate  of  ^  based  on  Y  such  that  jjCY-^i)  ■  £QQ  -  %  . 

A 


The  remark  shows  that  with  a  translation-invariant  location  estimate, 
it  is  sufficient  to  deal  with  the  case  ^  =  0. 

For  each  n  ■  1,2,  ...  let  ^  be  a  measurable  subset  of  n-dimensional 
Euclidean  space  and  consider  the  following  assumptions  for  the  sequence 


ASSUMPTION  (A1) 

sup  max  In-  |  0  as  n  •+•  »  . 

1  jn  .  in  1 

VMn  35 i<jsn 

ASSUMPTION  (A~)  For  any  6  >  0,  there  exists  integers  and  L  and 
partitions  Mn(l) ,  ...»  Mn(L)  of  ^  such  that  for  each  £=1,  ...,  L,  if 

d“«)  -  sup{Pj„-vin:  ^  .  %(*)> 

dij«>  *  ^  ‘ 

then  (1/n)  l±<. CdV^f)  -  d <  6  for  all  n^Nj. 

Two  other  assumptions  that  will  be  needed  are  as  follows. 

ASSUMPTION  (A^)  -j^ere  exists  constants  Nq  and  BQ  such  that 

”2Zi<jaij  s  Bo  for  a11  "  *  No- 

ASSUMPTION  (A^)  If  t'.(t)  denotes  the  cdf  of  the  difference  of  two  independent 
random  variables,  each  with  density  f,  then  G  has  a  density  g  *  G*  with 
g(t)  continuous  at  t  =  0  and  0  <  g(0)  <  ®  . 

Note  g(0)  =  /f^  ■  y  . 

To  study  the  asymptotic  properties  of  the  estimate  in  (1.1)  consider 


the  function 


Rn(t’VV  =  ^  -  2ty  . 

LBMA  2.1  Let  (Mn>  be  a  sequence  of  subsets  satisfying  assumption  (A^ 
and  let  assumptions  (A^)  and  (A^)  hold.  Then  for  any  e 
n  ■  1,2,...  and  any  real  number  t  >  0 

This  lemma  shows  that  vfr  Hn(t/*^\,rjn,Yn)/2t  converges  in  probability 
to  y  but  the  result  is  not  strong  enough  since  an  appropriate  t  must  be 
chosen  to  fit  the  data  and  the  nuisance  location  parameter  must  be 
estimated  when  is  not  zero.  The  following  theorem  shows  that  the 
convergence  in  the  lemma  is  uniform  and  this  will  allow  the  construction 
of  a  suitable  estimate  of  y. 

THEOREM  2,1  Let  {Mn>  be  a  sequence  of  subsets  satisfying  assumptions 
(A1)  and  (A2)  and  let  assumptions  (A^)  and  (A^)  hold.  Let  tQ  be  a  positive 
number.  Then 


sup0s  t  st0  IV-VV1 
V  Wn 


-*■  0  as  n  +  ®  . 


Suppose  that  ^  *  ^(Y^)  is  an  estimate  of  j ^  satisfying  the  following 


assumptions . 


ASSUMPTION  (A$)  For  all  n  =  1,2,.. 


LQL-tL.)  *  Hn'Xn)  '  Bn 
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ASSUMPTIQN  (Ag)  For  all  e  >  0  there  exists  a  sequence  of  subsets  (Mn) 
satisfying  (A^)  and  (A^)  and  an  integer  N  such  that 

P0^n  e  ^  >  1  -  e  for  all  n  >  N. 

THEOREM  2.2  Let  assumptions  (Ag) ,  (A4) ,  (Ag)  and  ( A &)  hold.  Let  tQ  be 
a  positive  number.  Then 

sup0stst  JV'-Kn-y  H  *°  • 

0  «n 

Let  the  n  x  1  residual  vector  be  and  the  empirical  cdf  of 

the  absolute  differences  in  ^  be  Hft)  =  J.  .a.  .<K|Z.  -  Z.|,t).  The  quantile 

i  1  J  ^ 

A  A 

of  order  a  for  vdll  be  denoted  by  t^.  Similarly  the  quantile  of  order  a 
for  H(t)  =  G(t)  -  G(-t),  t  >  0,  will  be  denoted  by  tQ  . 


LENMA  2.2  Let  as  suctions  (Ag) ,  (Ag)  and  (A&)  hold.  Then  for 
0  <  a  <  1, 

P 

-*■  t  as  n  +® 
a 

For  a  given  a  e  (0,1)  there  would  exist  tQ  so  that  za  s  tQ  with 
arbitrarily  high  probability  for  large  n.  Then  Remark  2.1  (iii) ,  Theorem 
2.2  and  Lemma  2.2  yield  the  main  result  that  follows. 


COROLLARY  2.1  Let  assumptions  (Ag) ,  (A4) ,  (Ag)  and  (Ag)  hold.  Let 
0  <  a  <  1  and 

Y  s  )/(2t J/R  )  . 

*  i  P 

Y  L  "*■  Y  as  n  -*■  «  . 

,un 

Hie  estimate  y  arises  from  yt.j^  #(3fn)  by  using  t  “  to  - 


Then 
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where  <  . . .  <  denote  the  M  =  (^)  ordered  absolute  differences 

*  Zj “Zi I  -  The  estimate  y*  of  (1.2)  arises  by  using  t  =  Sn  W^,  where  k=k^ 

converges  to  zero  at  an  appropriate  rate.  Specifically,  if  o  <  ax  <  a2 

and  kn  satisfies  a^  <  i^T  k^/M  <  ^ or  n»  then  is  bounded  away 

n 

from  zero  and  from  above  in  probability  and  Theorem  2.2  gives  the  following 
result. 

fflroUffiY  2.2  Let  assumptions  (Aj) ,  (A4),  (A$)  and  (A&)  hold.  Assume 
there  are  positive  constants  ap  &2  such  that  0  <  a.  <  £  k  /M  <  a  for 
all  n.  Let  n  2 

Y^kn/ZM^  . 

.  i  P 

Then  y*  ->  y  asn  +  ®  . 


.  APPLICATION  TO  LINEAR  MODELS 


In  this  section  it  is  shown  that  the  assumptions  of  the  previous 
section  are  satisfied  for  a  linear  model  under  mild  conditions  and  as 
a  result,  y  or  y*  can  be  used  as  a  consistent  estimate  of  y.  Thus  when 
using  rank  procedures,  a  consistent  estimate  of  y  is  available  without 
assuming  symmetry  for  the  error  distribution. 

Consider  the  model  (2.1)  with 


(3.D  ^  =  B01  +  X  A  /Jn  , 

where  1  is  an  n  x  1  vector  of  ones,  X  *  (x_)  is  an  n  x  p  design  matrix 
30(1  A  =  (Ai>  • • • »  Ap)  is  a  p  x  1  parameter  vector.  Let  aq  be  a  positive 
number  and  V  =  aq  for  k  *  1,  ...,  p}.  For  each  integer  n  =  1,2, 

let 

(3.2)  =  (0q1  +  X  A/ vft  :  a  c  P) 

ASSUMPTION  (B^)  For  each  k  =  1,  ...,  p, 

Iimlsi£n^xik"xk^>^  4  0  as  n  ->  ®  , 
where  is  the  average  of  the  kth  column  of  X 
ASSUMPTION  (B2) 

(l/n)X^.Xc  -*£  as  n-M»  , 

where  £  is  a  p  x  p  positive  definite  matrix  and  Xc  *  (I  -  (1/n)  J)^  is  the 
centered  design  matrix. 
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REMARK  3.1  Let  assumptions  (B1)  and  (B2)  hold.  Then  the  sequence  of 
subsets  Mn  in  (3.2)  satisfies  assumptions  (A^)  and  (A2) . 


Proof:  For  ^  t  Mn  write  njn-  r)in  -  xilPV'5  '  Then 

n“lSl<j5nl"jn-r'inl  5  2i0  2k.lmaxlsisntlxik  '  >•  if  *  e  V'  lhen 

(B^  implies  (Aj) . 

Next  let  6’  >0.  Then  there  is  a  partition  V,,  . . . ,  of  V  such  that 

for  all  1  =  1,  ...,  L,  A,  A’  e  implies  | Ak-A^.|  <6*  for  each  k  =  1 . p 

Then  Mn(£)  =  { BqI  +  X  a/^  :  A  e  t>£)  for  £  =  1 ,  . . . ,  L,  is  a  partition  of 

Mn  * 

Now 


(l/nJLjld^O-dJjli)]2 


•  t  s“?  .  iLitxjk- xik)(v ‘k5^' 


< 


(1/n 2)I 

i<j 


sup 

A,A'eP^ 


s  (P«,2/n2)lJ.1I1<j(xjk-  xik)2 

-  (PS'2/"2)!^!  n IJ.j  Cxik-  ^l2 

•  p«'2  trace  GSJ^/n)  . 

Since  5’  is  arbitrary  and  (B2)  holds,  it  follows  that  (A2)  holds.  □ 
Suppose  that  A  =  A(Y)  is  an  estimate  of  A.  The  intercept  parameter 
8q  has  no  effect  on  matters  here  since  it  can  be  absorbed  in  the  error 
variables  of  model  (2.1). 


ASSUMPTION  (Bj)  Let  A  satisfy  A(Y-X  &/Sn)  =  A(Y)  -  A  and  let  a  be 
bounded  in  probability  when  a  =  0  in  model  (3.1). 

This  assumption  is  satisfied  for  the  usual  rank  estimates  of  A 
proposed  by  Jaeckel  (1972),  Hettmansperger  and  McKean  (1977)  and  others. 

REMARK  3.2  Let  assunptions  (B^)  and  (B2)  hold  and  let  A  satisfy 
assumption  (B^) .  Then  assumption  (A^)  holds  for  in  the  column  space 

of  X  and  assumption  (Ag)  holds  with  Mn  given  by  (3. 2) . 

Proof:  Assumption  (A^)  follows  directly  from  (Bj) .  For  e  >  0,  the 
boundedness  in  probability  of  £  implies  that  a  sufficiently  large  aq  can 
be  chosen  in  defining  V  for  (3.2)  so  that  Pq(a  e  D)  >  1  -  e  for  large  n. 
Then  (A&)  follows  since  A  e  V  implies  6q1  +  X  e  Mn-  D 

These  remarks  show  that  (B^) ,  (B2)  and  (B^)  can  replace  (A,-)  and  (A^) 
in  Corollaries  2.1  or  2.2  in  the  case  of  linear  model.  Thus  under  typical 
assumptions,  y  or  y*  provides  a  consistent  estimate  of  y. 
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4.  LEHMANN'S  TWO- SAMPLE  ESTIMATE 


Lehmann  (1963)  proposed  a  consistent  estimate  of  y  for  the  two-sample 
problem  based  on  the  length  of  a  confidence  interval  for  a  shift  parameter. 
It  is  interesting  to  note  that  the  estimates  y  and  y*  are  computationally 
similar  to  Lehmann's  estimate  if  the  set  of  residuals  is  duplicated  to  form 
a  second  sample.  Lehmann's  work  showing  consistency  is  based  on  the 
assumption  of  two  independent  samples  and  would  not  apply  in  the  present 
context. 


To  describe  Lehmann's  estimate  let  two  samples  be  denoted  by 
*  •••»  Yjjj  and  X^,  ...,  X^.  Define  U(A)  =  4>(}C>Yj-A).  Define 

Ay  and  al  by  U(Ay) =  d^-Cl/2)  and  U(Ay)  =  mn  -  da+(l/2) ,  where  da  = 

(mn/2)  -  za^2(mn(m+n)/12)1^2  and  is  the  upper  a**1  quantile  of  the  standard 

normal  distribution.  Then  Lehmann's  estimate  of  y  is  • 

Consider  the  case  where  the  two  samples  are  identical,  that  is 

2 

m  =  n  and  Y.  =  X. .  Then  U(-A)  =  n  -  U(A) ,  A,  =  -A„  and  Lehmann’s  estimate 

11  L  U 

becomes 


(4.1)  (2za/2  //SH)/2Au  . 

Now  suppose  that  the  residuals  2^  of  section  2  are  duplicated  to  form 
two  "samples".  Use  equal  weights  a..  «=  l/(“).  Then  for  t  >  0, 

L 

H„(t)  -  i  -  cj)'1 

-  i  -  ("r'uctj  . 

Using  this  in  the  equation  defining  a^  yields  the  equation 
fin(V  '  ('1/n)  +  (2n/(n-l))za/2  /M.  . 
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Thus  for  large  n,  Ay  satisfies 
(4.2)  Hn(Ay)  =  2za/2  //6n 

and  with  (4.1)  it  follows  that  Lehmann's  estimate  is  approximately 
V  Ay)/2Ay  .  This  is  the  same  form  as  the  estimate  y  in  Corollary  2.1. 
Also  note  that  in  the  two -sample  model,  Lehmann  showed  /n  Ay  converges  in 
probability  to  a  constant  and  this  matches  the  behavior  of  t  in  y.  The 
probability  levels  a  used  in  each  estimate  are  not  related  in  any  useful 
way. 

In  Corollary  2.2,  suppose  kn  is  defined  by  kn/M  =  2za^2  •  Then 

(4.2)  shows  that  A„  corresponds  to  W,  and  the  estimate  y*  is  then 

Kn 

directly  in  the  form  of  Lehmann's  estimate  (4.1). 


-14- 


5.  APPENDIX 


Proof  of  LEMMA  2,1 
Write 

Rn(t’W  ’  Jj  aij  o'ij  - 2tY/ifi) 

where  Wi;j  =  (KjY.^-d^  | ,  t/*fi)  and  d„  =  rijn-nirj.  Using 

E0(Wij)  =  G(dij  +  t/y^  '  G(dij  •  t/y^ 

*  g(Cij)(2t/^rv) 

where  | j  _di j I  <  t /Jn  and  y  =  g(0),  it  follows  that 

’  2tJ.  aij(8«ii>-8«»)  • 

X<j  J 

Since  n  t  M,  Assumption  (A,)  implies  that  max  |d..|  ->■  0  and  then 
^  "  1  lsicjsn  1J 

max  |£. .  |  ->  0  as  n  -*■  The  continuity  of  g  at  zero  gives 
lsi<jsn 

|g(£^)  -  g(0)  |  <  e  for  all  lsi<jsn  for  any  t  >  0  if  n  is  sufficiently 
large.  Hence  for  any  e  >  0f|  Eg(Rn(t,r)n,Yn))|  <  2te  for  large  n  and  it 
follows  that  Eo^n^’^n’^n^  0  as  n  -*■  ®. 

Using  assumption  (A,)  and  |W..|  s  1,  a  standard  argument  shows  that 

O  lj 

var^CR^ft +  0  as  n  +  “>  .  With  the  mean  and  variance  tending  to 
zero  the  lemma  follows.  0 

Proof  of  THEOREM  2.1 


This  proof  is  somewhat  long  and  it  will  be  divided  into  pieces  with 
sane  lenmas.  Begin  by  choosing  any  e  >  0.  Choose  6  >  0  sufficiently  small 
so  that  4g(0)BQ*^  <  e/6,  where  BQ  is  specified  in  (Aj).  For  this  6  we  have 
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(take  Nj  >  Nfl  of  (Aj)),  the  partitions  1=1 . L  and  the 

extreme  differences  dV.(f)  and  df.  (l)  of  (A,).  Choose  <5 '  >  0  so  that 

lj  lj  L 

6g(0)6'  <  e/6.  Partition  the  interval  f  =  [0,tQ]  into  subintervals 
T1  ’  *  •  •  *  wi(ith  less  than  6 ' .  For  each  m  =  1 . M  let 

t!f  =  sup{t:t  €  T  }  and  t^  =  inf{t:t  €  T  }  . 
m  mm  m 


Then  for  each  m  *  1 ,  . . . ,  M 
(5.1) 

For  each  £=1,  ...,  L  and  m  *  1,  M,  define 


|t^  -  t^|  <  6’ 
1  m  m' 


Si j (i  ,m)  * 


i,  ifYj-Yj.  (dJjW 

or  Yj-Yj  <  (ijW  *  ,  d”  m*  J 

0,  otherwise 


for  l5i<j<n  and  Q.  m  *  ^fv  Y  a. -S.  ■  Gf.m)  . 

t,m  ij  vj 

LENMA  5.1  For  each  l  =  1 ,  ....  L  and  m  *  1 ,  . . . ,  M 


sup  -  yf.,0',301  S  Q*  *  w 

t  •  w  Cl 

m 

n,rj,«Mn(£) 

Proof:  Fix  l  and  m  and  write 


yt,n,Y)  -  yt»,n».x) 

-  vfi  l  a..(W..-W!.)  -  2y(t-t') 

where  W..  -  <K|Y.-Y.  -  d..|,  t/S),  W!.  -  <K|Y.-Y.  -  d!.  ( ,  t'f/H  )  , 
dy  *  hjn  -  and  d!  j  -  njn  -  .  Inspection  of  possible  cases  shows 


that  lwij‘w^jl  -  S^(£,jjO  for  all  lsi<jsn,  all  rj,  rj'c  Mn(£)  and  all 

t.t'  e.  T  .  Then  the  lemma  follows.  0 
m 

Using  Assumption  (A^) ,  there  is  6"  >  0  such  that  |d|  <  6"  implies 

g(d)  <  2  g(0).  Thai  from  assuiption  (Aj)  there  exists  N2  >  Nj  such  that 

Injn'^inl  <  5"/2  for  all  £  all  lsi<jsn  and  all  n  z  This  implies 

(5.2)  (-6"/2)  <  dE  (*)  s  dV.(£)  s  (6”/2) 

for  all  l  -  1,  ...,  L,  all  lsicjsn  and  all  n  s  N2. 

LE1WA  5.2  For  all  l  =  1 ,  . . . ,  L  and  m  =  1 ,  . . . ,  M, 

(a)  there  exists  an  integer  such  that 

E0(Q*,m)  5  (e/6)  +  4  g(0)6’  for  a11  n  *  N3’ 

(b)  var^CQ^  m)  •>  0  as  n  +  »  . 

Proof: 

E0^  =  *  l.  aij 

-  ^  .1.  aij < [GCdjj (A)  -  tj/yff  )  -  G(dEj(,)-t>  )) 

+  [G(dV.(t)  +  )-  G(dE.(f)+  tE/vfi  )]} 

-  &  _I_ai;j{g(Ci:j(£,ra))[d^(t)-  dEjW+  (tJJ  -  tE)/v£  ] 

-  g (Vy  (£ ,m) )  [dV. (t) -  dE.(f)+  (tJJ  -  tjp/yfi  ]} 

where  dE.  (1) -tjj/yfi  s  s  d^tf)-  tj/yff  and  d \.(l)  *  tE/^  s 

<  dVj  (£)  +  tjj/vfi-  for  all  lsi<jsn. 


Now  there  is  an  integer  Nj  >  ^  such  that  t^/ /n  <  6"/2  for  all  n  > 
This  along  with  (5.2)  implies 

-6"  <  Cij (£,m)  <  6" 

-6"  <  v>ij  (f,m)  <  6" 

for  all  ls;i<j<n  and  all  n  >  N^.  Then  by  the  choice  of  6”  it  follows  that 
s  *  4  sCOJ.I.a^TO-  .  ct“  -  t£)Mf  ] 

<  4  g(0)  &  l  a..[dj At)  -  dK (£)]  +  4  g(0)6* 
i<j  •LJ 

from  (5.1)  for  all  n  5  Nj,  Applying  the  Cauchy-Schwarz  inequality  and 
Assumptions  (.Z^) ,  (A^)  yields 

Efl(Q,  J  s  4  S(0)  [n2  I  af  (1/n)  I  (d "(f)  -  dj, <*) ) 2] 1/2+  4  g(0)6 

u  i<j  i<j  1J 

s4  g(0)  BQy^  +  4  g(0)  5' 
s  (e/6)  +  4  g(0)  6' 
from  the  choice  of  6  . 

The  second  part  of  this  lemma  follows  from  Assumption  (Aj)  and  the 

fact  that  \Stj  (£,m)  |  s  1  with  a  standard  variance  argument.  □ 

Choose  an  arbitrary  e’  >  0.  For  each  1*1,  . . . ,  L  and  m  *  1 . M 

/ 

choose  a  sequence  rj^  c  Mn(£)  ,  n  *  1,2,  ....  and  a  number  tffl  «  Tffl  . 

LEM4A  5.3  There  exists  an  integer  N4  such  that 

po(|wjS>vi  >  e/3)  s  e,/2ui 

for  all  l  *  1,  ...»  L,  all  m  ■  1,  ...,  M  and  all  n  *  N^. 

Proof:  Apply  Lenina  2.1  for  each  of  the  finite  number  of  pairs  (l ,m) .  Q 
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With  these  preliminaries  completed,  the  proof  of  the  theorem  will 
now  be  completed.  For  each  l  =  1 ,  . . . ,  L  and  m  =  1 ,  . . . ,  M  write 


teT 

m 

*  v  ?■?  -  w&vi  *  iw^oyi 2  e  > 

teT 

m 


P0«J>  * 

2^'  ♦  2  E  > 

w  - 

¥<V 

*  Eo «t,„>  +  IVcm 

4- V  l*o 

P0^,n,  ‘ 

w 

*  2 

2e/3) 

P0«<.m  ‘ 

E0 

2  *  V'W 

►nj#  V  1  * 

(9/c2)  Var0(Qi>in) 

♦  (c'/2LM)  , 

e'/LM 

where  Lemma  5.1  was  used  in  the  third  line.  Lemma  5.2(a)  and  the  choice 
of  6'  were  used  in  the  fifth  line,  Chebyshev's  inequality  and  Lemma  5.3 
were  used  in  the  seventh  line,  and  the  last  line  follows  from  Lemma  5.2(b), 
since  the  first  term  in  the  seventh  line  above  can  be  made  less  than  e */2LM 
for  all  l, m  for  n  sufficiently  large,  say  for  n  2  N.  >  N,.  Finally, 
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V  IV‘-VXn>lieJ 

VMn 


m 


VMnK> 


*  Ci4.i',/U' *  •* 


for  all  n  2  N5  and  this  completes  the  proof.  [] 

Proof  of  THEOREM  2.2 

With  assumption  (A,.)  and  Remark  2.1  it  is  sufficient  to  prove  the 
result  when  ^  =  0  .  Let  e  >  0.  Then 

P°(0SSS’ t  lWt’«n»Vl  >  €  ’ 

S  I'otO£SSt0lRhCt’£"’Xn)l  >  C’  *  V  *  P0(5n  1  Hn> 

S  Po<teSt  |Rn(t,fln,Xn)l  ”  c  5  *  P0<£„  *  > • 

VMn 

Then  Theorem  2.1  and  assumption  (A^)  imply  that  this  can  be  made  arbitrarily 
small  for  n  sufficiently  large  and  the  result  follows.  Q 

Proof  of  LENKA  2.2 

The  proof  is  basically  straightforward  and  brief  details  will  be  given 
here.  Assunption  (A&)  implies  that  it  is  sufficient  to  deal  with  the  case 


Kh  -  0.  Noting  that  Z^-l  =  (Y^ -Yi)  -  (Pjn-V'in)»  assunption  (Ag)  (only 

the  (A,)  part  is  needed)  implies  that  max  |y.  -y.  |  will  be  arbitrarily 
1  lsi<jsn  3n  in 

small  with  arbitrarily  high  probability  for  large  n  and  the  empirical  cdf 

of  Z.-2.  can  then  be  bounded  by  empirical  cdf's  of  Y. -Y. .  Assumption  (A_) 

J  1  J  1  3 

is  used  to  show  the  weighted  empirical  cdf  of  Yj -Y^  converges  in  probability. 
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